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Abstract. 

In p° photoproduction in ultra-peripheral heavy ion collisions, either ion can 
be the photon emitter or the target. The two possibilities are indistinguishable, 
and they should be able to interfere, reducing rho° production at low transverse 
momentum, pt < h/(b), where (b) is the median impact parameter. 

The two p° production points are separated by (6) m 18 — 46 fm, while the p° 
decay before travelling 1 fm. The two decay points are well separated in space- 
time, so the decays proceed independently and any interference must involve the 
final state n + n~ . This requires a non-local wave function. 
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Figure 1. (a) The two possibilities for exclusive p° production. These two 
diagrams may be accompanied by mutual Coulomb excitation, as in (b). Different 
time orderings are not shown in (b), but do not affect the interference. 



1. Introduction 

p° may be produced electromagnetically (via photoproduction) in distant interactions 
between ultra-relativistic heavy ions. In these ultra-peripheral collisions (UPCs) £Q, a 
photon from the electromagnetic field of one nucleus fluctuates to a quark-antiquark 
pair, which then scatters elastically from the other nucleus, emerging as a vector meson 

Elastic scattering is colorless, and can be described in terms of soft Pomeron 
exchange. For light mesons like the p°, the cross section rises slowly with the photon 
energy k. For small momentum transfers, the elastic scattering is coherent over the 
entire nuclear target, greatly increasing the cross section. For heavy mesons, the 
scattering cross section scales as the atomic number A 2 , while for lighter mesons 
like the p°, a Glauber calculation finds that shadowing reduces the dependence to 
approximately A 5 / 3 00]. 

The large photon flux and the coherent scattering lead to a large cross section for 
p° production. For 200 GeV per nucleon Au-Au collisions at RHIC, the cross section 
for p° production is 590 mb, or about 8% of the hadronic cross section The STAR 
collaboration has already measured p° production in Au-Au collisions at an energy of 
130 GeV per nucleon The cross sections agreed with the Glauber calculations. 

The elastic scattering is mediated by the strong force, and so has a short range 
0. The p production is localized to within 1 fm of the two ions. Figure la shows 
the two possibilities: either nucleus 1 emits a photon which scatters off nucleus 2, 
or vice versa. These two possibilities are indistinguishable, and are related by a 
parity transformation. Vector mesons have negative parity, so the two amplitudes 
subtract, with a transverse momentum (pr) dependent phase factor to account for 
the separation. The cross section is [7] 

I 1 2 

a = | A\ — A2 exp (ipT ■ b) | (1) 

where A\ and A2 are the amplitudes for p° production from the two directions. At 
mid-rapidity A\ = A2 and 

<7 = <7 (&)[(1-cos0t-&)] (2) 
where o~o(b) is the cross section without interference. The system acts as a 2-slit 
interferometer, with slit separation b Of course, b is unmeasurable, and the 
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observed pt spectrum is obtained by integrating Eq. (1) over all 6, subject to the 
constraint that there are no hadronic interactions (roughly equivalent to b > 2R A ). 
The cross section is suppressed for p? < h/(b), where (b) is the median impact 
parameter. Interestingly, for photo 

The p decay distance, 7/3ct < 1 fm is far less than (b) w 46 fm. The p° decay 
before the amplitudes from the two sources can overlap, making this an interesting 
test of quantum mechanics [H] • 

The coupling constant Za (where a « 1/137 is the electromagnetic coupling 
constant) is large, so multiple interactions between a single ion pair are common, the 
p° production may be accompanied by mutual Coulomb excitation, as in Fig. lb. 
Each ion emits a photon, which excites the other nucleus, usually to a giant dipolc 
resonance (GDR). The resonance decays by emitting one or more neutrons. These 
neutrons are a clear experimental signature for nuclear dissociation. 

The mutual Coulomb dissociation also 'tags' collisions with moderate impact 
parameters, 2R A < b <« 30 fm JUj. As long as k < jh/b, the photon flux drops as 
1/b 2 . For a single photon reaction like p° production, the production probability P(b) 
follows the same scaling. For n-photon interactions, P(b) ~ l/b 2n ^Tj so the median 
impact parameter is much smaller and interference is visible at larger px- 

Except for the shared impact parameter, the photon exchange reactions are 
independent The cross section for p° production with mutual Coulomb 

dissociation is 



where P p (b) and P2EXc(b) are the probabilities for p° production and mutual 
excitation respectively. Here, Phad(b) is the probability of having a hadronic 
interaction; the last term has an effect similar to setting a minimum impact parameter 
b mm = 2R A . For gold ions, P p (2R A ) « 1% and P 2 exc(2R a ) « 30%, so the 
probability of multiple interactions is substantial. Through the smaller bs, the photon 
spectra in multi-photon interactions are harder than for single photon reactions. 

2. The STAR Detector 

The Solenoidal Tracker at RHIC (STAR) studies heavy-ion and polarized proton 
collisions at the RHIC collider Brookhaven National Laboratory. STAR is optimized 
for high-multiplicity final states, but it is also effective for 2-4 charged particle final 
states from UPCs. STAR collected UPC data in gold-gold collisions at energies of 
130 GeV per nucleon (in 2000) and 200 GeV per nucleon (in 2001). STAR has also 
studied p° photoproduction in 200 GeV per nucleon deuteron on gold collisions |T2|. 
The p° interference analysis uses the 2001 gold-gold collisions dataset. 

Charged particles are reconstructed in a 4.2 meter long, 4 meter diameter 
time projection chamber (TPC) ^H]- This data was taken with the TPC in a 0.5 
T solenoidal magnetic field. The tt^ reconstruction efficiency is high for tracks 
with transverse momentum px > 100 MeV/c and pseudorapidity < 1.15. The 
track position and specific energy loss, dE/dx were measured at 45 points for high 
momentum charged particles with \r}\ < 1. The TPC is surrounded by 240 scintillator 
slats covering \rj\ < 1, comprising the central trigger barrel (CTB). Two zero degree 
calorimeters (ZDCs) are located 18 m upstream and downstream of the interaction 
point. These calorimeters are sensitive to neutrons from mutual Coulomb dissociation 




(3) 
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Figure 2. The (a) rapidity and (b) tttt invariant mass spectra of p° using data 
from the minimum bias trigger. The mass is fit to to a Breit-Wigner function for 
the p°, plus a Soding term due to interference with a direct tttt component. The 
combinatoric hadronic background is also shown. 



Two different triggers jTH] were used to study interference. The 'topology' trigger 
selected events with 2 roughly back-to-back tracks in the CTB. It divided the CTB 
into 4 quadrants: north, south, top and bottom, and required hits in the north and 
south quadrants. The top and bottom regions were used as vetoes to reject cosmic 
rays. 

The minimum bias trigger selected events with one or more neutrons in each 
ZDC. The ZDC signals were required to occur within 1 nsec, restricting these events 
to the central 30 cm of the TPC. Data from both triggers was processed identically, 
except that events from the CTB based trigger were distributed more broadly along 
the TPC axis, and consequently, were accepted in a broader range. 

3. p° Production 

UPC p° production has a distinctive signature - two oppositely charged tracks with 
small net pr- We select events with exactly 2 primary tracks that form a vertex. The 
vertexing procedure considers all of the tracks in an event, and rejects tracks that are 
inconsistent with coming from a single vertex. We allow a few non-primary background 
tracks in the TPC. The topology and minimum bias data were treated identically, 
except that allowance was made for the different distribution of the accepted event 
production points. This study used about 1.5 million minimum bias and about 1.7 
million topology triggers. This analysis will study the pt spectra of produced p° . 

Figure 2a shows the rapidity distribution of the coherent p° (pairs with px < 150 
MeV/c). The points (data) are in excellent agreement with a calculation based on the 
soft Pomeron model and our detector simulation. 

Figure 2b shows the tttt invariant mass, M Ilr for the sample. The data is fit to a 
relativistic Breit-Wigner for the p°, plus a Soding term to account for the interference 
with direct ir + ir~ production [lfi| . The interference shifts the peak of the distribution 
to lower M^tt- The direct tttt to p ratio agrees with the STAR 130 GeV analysis [1] 
and with studies by the ZEUS collaboration in jp interactions |17j . However, with 
the nuclear coherence, the STAR data is at smaller |f | than the ZEUS analysis, so the 
two results may not be directly comparable |18j . 
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Figure 3. Raw (uncorrected) t± spectrum for p° sample for 0.1 < \y\ < 0.5 for 
the topology data. The points (red) are the data. The dashed histogram (blue) is 
a simulation that assumes that there is interference, while the dot-dash histogram 
(black) is based on a calculation without interference. The solid black histogram 
with very few counts is the like-sign background. 

4. p° Interferometry 

Stringent event selection criteria were used to select a clean, low-background sample 
for the p° interference study. Events were required to have exactly two tracks with a 
vertex within 50 (100) cm longitudinally of the center of the TPC for the minimum 
bias (topology trigger) sample. p° candidates with \y\ < 0.1 were eliminated to 
avoid possible contamination from cosmic rays. The pairs were required to have 
550 MeV < Mtttt < 920 MeV. The minimum mass cut removed most of the hadronic 
backgrounds, which are concentrated at low M^. These cuts reduced the like-sign pair 
background to a few percent. With these cuts, the background from misidentified two- 
photon reactions like 77 — + e + e _ should be very small. The sample still includes direct 
pions, which should have the same spin/parity and quantum mechanical behavior as 
the pion pairs from p decay. We do not distinguish between the two sources. 

The magnitude of the interference depends on the ratio of the amplitudes for p 
production on the two nuclei. Away from y = 0, the photon energies for the two photon 
directions differ, ki,2 = My/2exp(±y/2); the amplitudes differ and the interference is 
sub-maximal. Although it is not expected in the soft-Pomeron model, the the photon 
energy difference could introduce a small p° production phase difference, which could 
affect the interference j2J- This analysis focuses on the region near mid-rapidity where 
any phase difference should be small. A Monte Carlo is used to find the interference 
for different rapidity ranges • 

We study the pr spectra using the variable t± = p\. At RHIC energies, the 
longitudinal component of the 4- momentum transfer is small, so f « t±. Without 
interference, the spectrum dN/dt oc exp (— bt) [71 ITT?]. 

Figure 3 compares the uncorrected minimum bias data for 0.1 < \rj\ < 0.5 with 
two simulations, with and without interference. Both simulations include the detector 
response. The data has a significant downturn for t < 0.001 GeV 2 , consistent with 
the (b) = 18 fm expected for a p° accompanied by mutual excitation JTj. This data 
matches calculation with interference, but not the no-interference result. 

The efficiency corrected data are shown in Fig. 4. Minimum bias and topology 
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Figure 4. Efficiency corrected t± spectrum for p° from (a) mutual dissociation 
with 0.1 < \y\ < 0.5, (b) mutual dissociation with 0.5 < \y\ < 1.0, (c) topology 
trigger with 0.1 < \y\ < 0.5 and (d) topology trigger with 0.5 < \y\ < 1.0. The 
points arc the data, while the solid line is the result of the fit to Eq. 4. 



data are shown separately, each with two rapidity bins: 0.1 < \y\ < 0.5 and 
0.5 < \y\ < 1.0. The efficiency is independent of pr- However pr smearing (resolution) 
does affect the spectrum slightly; the efficiency correction removes this. The p° pr 
resolution is about 9 MeV/c, while the 1st t bin covers to (15 MeV/c) 2 . Interference 
depletes the first few bins, but feed down from the higher t bins partially repopulates 
them. 

The data is fit to the 3-parameter form: 
rIN 

^— = a exp(-fa) [1 + c(R(t) - 1] (4) 

where R(t) = Int(t)/Noint(t) is the ratio of the Monte Carlo t-spectra with and 
without interference. Here, a is the overall normalization, the slope b w R 2 Al and c is 
the degree of spectral modification; c = corresponds to no interference while c = 1 
is the calculated interference. This functional form separates the interference (c) from 
the nuclear form factor (b). 

Table 1 gives the results of the fits. At small-rapidities the amplitudes from the 
two directions are similar and the interference reduces the cross section at p T = by 
more than at larger rapidities. In the minimum bias data, the interference extends to 
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Trigger |Rapidity| b (GeV~ 2 ) c x 2 /DOF 

Minimum Bias 0.1 - 0.5 301 ± 14 1.01 ±0.08 50/47 

Minimum Bias 0.5 - 1.0 304 ± 15 0.78 ± 0.13 73/47 

Topology 0.1 - 0.5 361 ± 10 0.71 ±0.16 81/47 

Topology 0.5 - 1.0 368 ± 12 1.22 ±0.21 50/47 



Table 1. Results of the fits to the two datasets for the two rapidity bins. 



higher p? than the topology data because the former has a smaller (b) . 

The 4 c values are consistent within errors; the weighted average is c = 0.93±0.06. 
The b values for the minimum bias and exclusive p° data differ by 20%: 364±7 GeV -2 
for the exclusive p versus 303 ± 10 GeV -2 for the Coulomb breakup events. 

The different b values may be attributed to the different impact parameter 
distributions caused by the nuclear breakup tagging in the minimum bias data. The 
photon flux depends on the impact parameter as l/b 2 . When b s=s fewi?^, p are more 
likely to be produced on the side of the target near the photon emitter than on the 
far side. p° production is concentrated on the near side, leading to a smaller effective 
production volume and the smaller b. This near-side skewing affects the interference 
slightly, but is not included in current calculations. 

Systematic errors come from a variety of sources. We have studied the effect of 
detector distortions (primarily px smearing) by turning off the detector simulation and 
comparing raw simulations with reconstructed data. This lowered c by 18%; if the 
detector simulation is 75% correct (a very conservative assumption), then detector 
effects are a less than 5% systematic uncertainty. Minor systematic uncertainties 
come from backgrounds and the fitting procedure. We have compared the data and 
simulations for a variety of kinematic and detector-based variables, and found good 
agreement. We estimate a preliminary experimental systematic uncertainty of 8%. 

This analysis depends on the shape of R(t), which is calculated following Ref. 
There are some simplifications in Ref. jxj. The calculation averages the photon flux 
over the surface of the nucleus, rather than incorporating the proper 1/b 2 weighting. 
This may partly explain the relatively poor fit for the large- \y\ Coulomb breakup 
fit. We estimate that the effective impact parameter should be within 10% of the 
actual ion-ion separation for the coulomb breakup data, and 3% for the exclusive 
p. The uncertainties in the calculations should be at most a 15% effect. With this, 
the interference is 93 ± 6(stat.) ± 8(syst.) ± 15(theory)% (STAR preliminary) of that 
expected. 

The p decays rapidly, with ~/f3cT <C (b) and the two p decay points are 
well separated in space-time. In the usual space-time picture, the decays occur 
independently, and any interference must involve the final state 7r + 7r~. Interference 
must involve identical final states from the two sources. However, given the large 
available phase space for the decays, this is very unlikely for independent decays. 

One interpretation of this result is given in Ref. In it, the interference 

occurs because the post-decay wave function includes amplitudes for all possible final 
states, then the amplitudes for identical states subtract, and the interference is visible. 
Because of the two sources, the tt + tt~ wave function is non-factorizable, and thus 
exhibits the Einstein-Podolsky- Rosen paradox [201 - We nn d that the decoherence, 
1 — c, due to environmental or other factors is less than 43% at the 90% confidence 
level. 
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5. Conclusions 

The STAR collaboration has studied exclusive p° photoproduction and p° 
photoproduction accompanied by nuclear excitation. The cross sections and rapidity 
distributions match the predictions of the soft Pomeron model. 

Production at low— pt is suppressed, as expected if p° production at the two ions 
interferes destructively. The suppression is 93 ± 6(stat.) ± 8(syst.) ± 15(theory)% of 
that expected. The decoherence is less than 43% at a 90% confidence level. 

We thank Tony Baltz and Joakim Nystrand for help with the theoretical 
calculations. This work was supported by the U.S. DOE under contract number 
DE-AC-03-76SF00098. 
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